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ABSTRACT 

We present the discoveries of three faint companions to young stars in the Scorpius-Centaurus region, 
imaged with the NICI instrument on Gemini South. We have confirmed all three companions through 
common proper motion tests. Follow-up spectroscopy has confirmed two of them, HIP 65423 B and 
HIP 65517 B, to be brown dwarfs, while the third, HIP 72099 B, is more likely a very low-mass star just 
above the hydrogen burning limit. The detection of wide companions in the mass range of ^40-100 Mj up 
complements previous work in the same region, reporting detections of similarly wide companions with 
lower masses, in the range of ^10-30Mj up . Such low masses near the deuterium burning limit have 
raised the question of whether those objects formed like planets or stars. The existence of intermediate 
objects as reported here could represent a bridge between lower-mass companions and stellar companions, 
but in any case demonstrate that mass alone may not provide a clear-cut distinction for the formation of 
low-mass companions to stars. 

Subject headings: brown dwarfs — planetary systems — binaries: general 



Introduction 



Although a number of brown dwa r fs are known to 
orbit stars (e.g. Nakajima et al.|[l995l iBurgasser et al 
20001). most detected brown dwarfs are single (e.g 



(UCL) and Lower Centaurus Crux (LCC). Age esti- 
mates of t he various sub-region s have varied in th e 



literature dde Zeeuw et al.l 1 19991 : IPecaut et al.l 1201 2|) 



Kirkpatrick et al. 2012) o r acc ompanying other brown 
dwarfs (e.g. iBouy et al.l 12008 ). In general, brown 
dwarf companions to stars ar e rare, particularly at 
small orbits of a few AU (e.g. iGrether & Lineweaver 
2006). However, indications exist that they may be 
at least slightl y more common in wider configu rations 
(e.g. lLafreniere et al.ll201 ltllreland et alj|201 lb . 



but recent direct comparisons to other young groups 
have shown that the whole region is younger than the 
[3 Pic moving group (12 Myr), placing the UCL and 
LCC sub-groups at ages of approxi mately 10 Myr, an d 
USco probably somewhat younger (ISong et al.ll2012l) . 
Previous detections of wide (^200-60 AU) low-mass 



subst e llar companions in USco (e.g. lLafreniere et al 



2008; lire land et al.l 1201 II) in shallow adaptive optics 



The Sco-Cen region (Ide Zeeuw et al .1 119991) is a 
young stellar association consisting of the sub-regions 
Upper Scorpius (USco), Upper Centaurus Lupus 
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(AO) surveys motivated us to perform an extended sur- 
vey in the whole Sco-Cen region using NICI at Gemini 
South. In the course of this study, we detected several 
new companion candidates that are being followed 
up for astrometric and spectroscopic confirmation. In 
this Letter, we report the discovery of three brown 
dwarf candidates, all of which have been confirmed 
though common proper motion and two of which have 
also been spectroscopically confirmed to be substellar. 
Two of the host stars, HIP 65423 and HIP 65517, are 
members of LCC, and HIP 72099 is a me mber of UCL 
dde Zeeuw et al]ll999t IChen etal.ll201lb HIP 72099 
has a known debris disk dChen et al.ll201lh . 
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2. Observations and Data Reduction 



2.1. Imaging Observations 

Images were taken with NICI at Gemini South as 
part of a program to search for planet and brown dwarf 
companions to stars in the Sco-Cen region. The dual 
band imaging mode was used, with the broad K s filter 
(2.15 /im with 14.5% bandwidth) in the red channel 
and the narrow H2 filter (2.12 /im with 1.23% band- 
width) in the blue channel. Each target was observed 
with a short exposure of 0.38 s integration time and 10 
coadds, and a longer exposure of 80 s in a single in- 
tegration. This setup was used to allow for a wide dy- 
namic range with sensitivity to bright and very close-in 
companions in the short exposure and narrow band, as 
well as faint and wide companions in the long expo- 
sure and broad band. The targets presented here were 
each first observed in April or June 201 1 and then re- 
observed on both 2012 Feb 29 and 2012 Apr 2. 

The data were reduced with a custom IDL pipeline. 
Flat fielding and dark subtraction was applied to all 
individual frames. While no dithering was used, ob- 
servations of different targets in the survey were taken 
contiguously and on different random places on the de- 
tector, thus forming an equivalent of the often used 'jit- 
tering' strategy where the science frames of other tar- 
gets were used as references for a given target frame, 
for sky subtraction and identification of bad pixels. All 
frames were then distortion corrected and re-oriented 
to a common orientation with North up and East to 
the leffl Ensemble astrometry showed that a system- 
atic 0.3° offset was present during one specific NICI 
mounting in June 201 1, so since one epoch of imaging 
of both HIP 65517 and HIP 72099 was taken during 
that period, those images have been corrected to ac- 
count for this offset. 

Reduced images for the three systems are shown 
in Fig. 03 The limiting magnitudes range from 19.4 
to 20.2 mag in the 80 s broad-band images, and no 
tertiary candidates were found in either system. 

2.2. Spectroscopic Observations 

Spectra for the three targets were acquired using 
NACO at the VET in May 2012 for HIP 65423 and 
HIP 65517 and July 2012 for HIP 72099. We used a 



1 By default, the NICI images are rotated with North pointing down- 
wards, and since the NICI port varied between up-looking and side- 
looking from 201 1 to 2012, the orientations varied between a right- 
handed and a left-handed configuration. 



m 




Fig. 1. — Discovery images of the three sys- 
tems. The projected separations range from 0.350"for 
HIP 65517 B to 1.835"for HIP 65423 B. North is up 
and East is to the left. Note that apparent point sources 
closer to the stars are due to known ghosts. 
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slit of 172 mas width and 40" length in the H + K set- 
ting giving spectra over a wavelength range of ~ 1 .33— 
2.53^m simultaneously at a spectral resolution of 550. 
Both the primary and secondary were included in the 
slit for all exposures. The observations were executed 
using an ABBA nodding pattern with a nod throw of 
10". One AB pair was taken with a short integration 
time of 1 s with 12 coadds for getting the primary star 
in a fully linear count regime, followed by two full 
ABBA cycles with longer integrations of 5 s with 20 
coadds per frame, for a total of 800 s on each target. 

Data reduction was performed using a custom IDL 
pipeline based on a previously develop ed code for 
NACO spectroscopy djanson et al.l l201oh . The data 
were fiat fielded and dark subtracted, and each contigu- 
ous AB or BA set was pairwise subtracted to remove 
the sky. For each individual frame, a trace of the spec- 
trum of the primary star was acquired through gaussian 
centroiding on each individual pixel row across the 
wavelength direction. This worked well, even though 
mild saturation occured in the iJ-band range for the 
long exposures. The traces were then re-centered us- 
ing spline interpolation to form a straight vertical trace 
centered on a common pixel column for all frames, af- 
ter which the short exposures were coadded to form a 
single non-saturated primary star reference frame, and 
likewise the long frames were coadded to form a single 
frame with maximal sensitivity to the companion. 

Three ID spectra are extracted from the spectral 
trace image, each in an aperture of 3 pixels: One 
centered on the primary (this is done on the basis of 
the short-exposure coadded frame, since the primary 
star is partly saturated otherwise), one centered on the 
companion, and another centered on the exact oppo- 
site side of the primary, to represent the PSF at the 
location of the companion. Wavelength solutions for 
these spectra are acquired using both sky lines and in- 
trinsic lines in the target spectra. The PSF reference 
spectrum is then subtracted from the companion spec- 
trum to form a PSF-subtracted companion spectrum. 
This in turn is divided by the primary spectrum, ap- 
propriately scaled to account for the difference in in- 
tegration time between the short and long exposures, 
which gives a telluric-free contrast spectrum. By mul- 
tiplying the contrast spectrum with a model spectrum 
for the primary, this yields a companion spectrum that 
is telluric-free and subject only to the uncertainties of 
the primary star properties. In our case, we simply 
employ a blackbody spectrum at the stellar tempera- 
ture as a model for the stellar flux, since the intrinsic 



stellar lines are weak in this wavelength range and do 
not affect the purely continuum-based analysis that we 
perform in Sect. 13.21 The adopted temperatures are 
5830 K for HIP 65423 (G3V), 5860 K for HIP 65517 
(G2V), and 6360 K for HIP 72099 (F6V), using SIM- 
BAD spectral types. 

3. Analysis 

3.1. Astrometric Analysis 

The main limiting issue for accurate astrometry 
for the NICI observations is that the minimum inte- 
gration time is 0.38 s, which for the bright stars ob- 
served in our program means inevitable saturation. 
Furthermore, the NICI PSF is strongly asymmetric 
and the non-linearity response of the detector is com- 
plex. This means that neither Moffat fitting of the PSF 
halo nor centroiding on the saturated core, which often 
work well for s i milar instruments like HiCIAO (e.g. 
Thal mann etall l201lh or NACO (e.g. Ijanson et al 



2008), can be satisfactorily applied here. Furthermore, 
the PSF spider arms, which have been succ essfully 
used for e.g. NIRI (e.g. lLafreniere et al . 2007) are not 
concentric for NICI. Yet, it is fairly straightforward 
to determine where the actual center is by eye. We 
have tried all these techniques on ~50 target stars with 
suspected background stars in the program, and found 
that manual centering gave the smallest scatter in mea- 
sured versus expected background motion, meaning it 
is likely the most accurate method. Hence, we use that 
for HIP 65423, HIP 65517 and HIP 72099 as well. 
As we will see below, this indeed leads to a good as- 
trometric quality. For astrometry of the companion 
candidates, we used Gaussian centroiding. 

The astrometry of each candidate relative to the ex- 
pected background motion is shown in Fig. [2] The 
scatter is consistent in the x- and y-positions for all tar- 
gets, and has an average of 5 mas, which we adopt as 
the astrometric uncertainty. It can be seen that each of 
the candidate companions is well consistent with com- 
mon proper motion. For HIP 65423, the significance is 
9.4cr, for HIP 65517 it is 1 l.Ocr, and for HIP 72099 it 
is 4.0(7. We thus consider all of them to be confirmed 
companions. 

3.2. Mass estimation 

Our primary objective for the spectroscopic anal- 
ysis is to generate mass estimates for the compan- 
ions. We do this in two independent ways: First, we 
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Fig. 2. — Astrometry for each system. Symbols with error bars are the measured positions. The solid line shows the 
predicted motion for a static background star, with circles4narking the expected location at each respective date. The 
background hypothesis can be excluded in each case, confirming common proper motion. 



calculate H- and if-band brightnesses for each com- 
panion. This is easily done since we have a direct 
contrast spectrum between the primary and compan- 
ion. Hence, we can calculate the contrast within a 
given band as the mean of the contrast spectrum within 
that band (weighted appropriately by the spectral re- 
sponse curve for H and K, respectively), and then 
add the measured H- and if -band magnitudes for the 
star to get the corresponding values for the compan- 
ion. Note that extinc tion is negligible for all three 
stars ( Chen et al1l201 1 ). Following this procedure, and 
converting to absolute magnitudes using Hipparcos 
( Perryman et al.l[l997 : van Leeuwenll2007h parallaxes, 
gives Aijj = 7.46 mag and Mk; = 7.04 mag for 
HIP 65423 B, Ai H = 7.63 mag and M K = 7.32 mag 
for HIP 65517 B, and M H = 6.90 mag and M K = 
6.58 mag for HIP 72099 B. Using th e BCAH evo- 
lutionary ^jn^dd^jBaraffeeLalJ LL998|) for an age of 
10 Myr dSongetal.ll2012h gives a mass of -65 M jup 



consistently in both H and K for HIP 65423 B. 
For HIP 65517 B, the mass is — 55Aij up and for 
HIP 72099 B the mass is ~95M jup consistently in 
both filters. 

For the second method, we compare the spectral 
continuum with library standard stars in order to deter- 
mine a spectral type (SpT) and temperature. The stan- 
dard stars are taken from the IRTF spectral library fo r 
M-type stars (ICushing et al .1l2005HRavner et al.ll2009l 



From visual inspection, it is clear that the compan- 
ions are in the SpT range of ~M6, so for a formal 
test we use the standard stars Gl 268AB (M4.5V), 
Gl 51 (M5.0V), GJ 1111 (M6.0V), Gl 406 (M6.5V) 
and Gl 644C (M7.0V) to encapsulate the feasible phys- 
ical range (there is no M5.5V standard in the library). 
These standard stars are probably older than our de- 
tected brown dwarfs, hence the surface gravities dif- 
fer modestly, but the two classes of objects are ob- 
servationally much clos er to each other th an to gi- 
ants or sub-giants (e.g. ISlesnick et al.l 120041) . and we 
have verified that such standard stars provide a visu- 
ally worse match to our companions than the lumi- 
nosity class V objects. We fit each standard star to 
each target spectrum and test the quality of the fit us- 
ing a squared residual (effectively \ 2 ) metric. We fit 
the H- and if -band ranges separately (1.49-1.83 /im 
and 1.97-2.33 /im, respectively), and re-normalize the 
target and standard to have the same mean flux value 
in the fitted range. In all cases (for each target and 
each wavelength range), the squared residuals as func- 
tion of spectral type appear to be continuous func- 



tions with well-defined minima, hence we fit a second- 
order polynomial and adopt the spectral type at the 
minimum of the fitted function as the best fit. For 
HIP 65423 B, this gives M6.03 in H and M6.31 in 
K, for HIP 65517 B, it gives M5.30 in H and M5.45 
in K, and for HIP 72099 B it gives M5.85 in H and 
M6.50 in K, all nicely consistent and in accordance 
with visual inspection. For translating spectral types 
into t emperatures, we use two different empirical re la- 
tions dGolimowski et alj|2004t ISlesnick et al.ll2004 to 



represent the uncertainty in this translation. The differ- 
ence between the relations is — 100 K at these spectral 
types. For s ome of the cases, we hav e to slightly ex- 
trapolate the Golimowsk i et al. (2004 ) relation, since it 
is calibrated for >M6.0. As best-fit values we consider 
the mean of the spectral types derived from H and K, 
and the mean of the resulting temperatures from the 
two relations. This gives 2770 K for both HIP 65423 B 
and HIP 72099 B, and 2960 K for HIP 65517 B. Using 
the BCAH model, we thus find a mass of ~45Aij up 
for HIP 65423 B, ~70M jup for HIP 65517 B, and 
~45M jup for HIP 72099 B. The spectra are shown in 
Fig. a 

Each of the methods used above have their respec- 
tive uncertainties due to significant uncertainties in the 
underlying parameters. For the brightness-based esti- 
mate, the dominant uncertainty is in the distance es- 
timation - due to the small parallax signature of these 
stars at their > 100 pc distances, the error is ~ 10-20%, 
leading to 0.3-0.5 mag uncertainties in the distance 
modulus. Our final upper and lower limits on compan- 
ion mass from this method also include photometric 
uncertainties, as well as the uncertainty in age, where 
we have adopted an age range of 8-12 Myr, given that 
the age of Sco-Cen has been deter mined as being be - 
tween that of TW Hya and /3 Pic dSong et al.l l2Q12h . 
The values are provided in table Q] 

The temperature-based estimate has the advantage 
of being independent of distance, but has its own un- 
certainties. For our error limits quoted in table [Tj 
we have included the scatter in spectral types de- 
rived from the different bands, the scatter in spec- 
tral type/temperature relation, and the uncertainty 
in age. One possible issue that may be relevant to 
keep in mind for the temperature-based estimation 
is the fact tha t in the eclipsing brown dwarf binary 
2M0535-05 dStassun et al.ll2006l) . the primary (more 
massive) component is colder than the secondary one. 
The reason for this effect is still under discussion, 
where arguments for and against physical mechanisms 
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Fig. 3.— Spectroscopy for HIP 65423 B, HIP 655 17 B, and HIP 72099 B. The target spectrum in black is wedged 
between the two best fitting spectral templates from lCushing et al.l (12005b and lRavner et al. I d2009h . 
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such as chromospheric activity have been pre sented 
dStassun et all 120121: iMohanty & Stassunll2012h . Our 
detected brown dwarfs are older than the ~1 Myr 
2M0535-05 system, so if the effect is age-related (e.g. 
due to activity or accretion), we may expect that the 
effect should be smaller in our case. Whatever the 
mechanism, the change in temperature is compensated 
for by an inflated radius in 2M0535-05 A, which keeps 
the luminosity largely constant, hence our magnitude- 
based estimates should be largely insensitive to this 
effect. In any case, the estimated masses for all three 
companions are consistent between the two methods, 
within the error bars. 

4. Discussion 

There is a well-known deficit of companions to 
stars with mass ratios in the range of ~ 1—10% at 
small semi-major axes (a few AU or less); the so- 
called brown dwarf desert (e.g. lGrether & Lineweaver 
2006). This desert is a natural consequence of the dif- 



ficulty for stellar systems to form at small mass ratios 
(e.g. [Bate 2009), and the simultaneous difficulty for 
planetary systems to form at large mass ratios (e.g. 
Mordasin i et al. 20091) . It is still unclear to which ex- 
tent the brown dwarf desert extends also to large semi- 
major axes (tens or hundreds of AU). While there are 
several indications that brown dwarfs are rare at these 



separ a tions (e.g . Kasper et al.l l2007t iNielsen & Close 
2010; iJanson et al.l 120121) some evidently do exist. 



Our detected brown dwarfs have best-fit mass ratios of 
4-5%, right in the middle of the desert. Other known 
systems with similar mass ratios and semi-major 



axes include e.g. GQ Lup (Neuhauser et al. 2005 



Jans on et alJl20 06: McElwain et al. | l2007l) and GJ 758 
dThalmann et al.l2009tlCurrie et al.l2010l;|janson et al 



2011) at 3—4%. While small compared to most bi- 
nary systems, these are nonetheless an order of magni- 
tude higher than the mass ratios of dir ectly imaged 
plane t ary sy stems, suc h as HR 8799 dMarois et al.l 
l2008i 1201 oh or /3 Pic dLagrange et al.l 120091 1201 Oh . 
Another interesting comparison is provided by the 
wide systems with companions of lower mass than 
our brown dwarfs that have been recently reported in 



the same yo ung region: 1 RXS J1609 dLafreniere et al.. 
20081 1201 Ol) HIP 78530 dLafreniere etal.l l201ll) and 
GSC 0621 dlrelandetailbOllh . with mass ratios of 



1-2%. Such low mass ratios by themselves may imply 
that these systems could have formed in a planet-like 
manner, perhaps through formation at a small semi- 
major axis and subsequent scattering outward through 



gravitation al interactions wi th other companions in 
the system dVeras et alj|20()3, e.g .). However, the fact 
that systems with intermediate mass ratios exist, such 
as those discovered here, at some non-negligible rate 
(perhaps ~3% if we naively divide the objects con- 
firmed so far by the number surveyed) complicates 
such an analysis. Although mass ratio remains highly 
useful to distinguish between formation scenarios, it 
is not as clear-cut for wide companions as it is closer 
to the star. Capture of free-fl oating objects in young 
stellar cluster environments dParker & Ouana 12012 : 
iPerets & Kouwenhovenl |2012|) is a mechanism that 
could potentially explain an extensive range of mass 
ratios at wide separation. A reservoir of such objects 
to be captured evidently exists, down to masses be- 
low the deuterium burning limit, though possibly not 
in sufficient num bers (e.g. iPena Ramirez et alj 12012 : 
Scholz et alJEoLi . 
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Table 1 : Observed and derived parameters for the dis- 
covered systems. 



Property 


HIP 65423 A 


HIP 65423 B 


HIP 65517 A 


HIP 655 17 B 


HIP 72099 A 


HIP 72099 B 


SpT 

T eff a (K) 

Dist. (pc) 
App. H (mag) b 
App. K (mag) b 
Mass c 

Ang. sep. (") 
PA(°) 

Proj. sep. (AU) 


G3 
5830 
124.4±26.6 
8.18±0.06 
8.08±0.02 
1.3 M Q 


M6 

27701™ 

12.93±0.06 
12.51±0.02 
65±|g/46llg M jup 
1.835±0.005 
247.4±0.2 
228±49 


G2 
5860 
110.6±18.6 
8.18±0.03 
8.08±0.03 
1.3 M Q 


M5.5 

2860±U0 

12.85±0.03 
12.54±0.03 

55tl° /70tf 5 M jup 
0.350±0.005 
321.7±1.4 
39±7 


F6 
6360 
157.7±40.6 
8.43±0.02 
8.40±0.03 
1.4 Mq 


M6 
2770jlgo 

12.89±0.02 
12.57±0.03 
95l^/45l 1 1 5 5 M jup 
0.667±0.005 
34.4±0.4 
107±27 



"Inferred from spectral types using ISherrvetalJ j2004) for the stars andjsiesnick et ai H2004h for the companions. 
6 Stellar magnitudes from 2MASS i Skrutskie et al. 200^). 

c Stellar masses inferred using D'Antona & Mazzitelli ( 1997), brown dwarf masses using Baraffe et al. ( 1998). For the brown dwarfs, the first value 
is based on intrinsic brightness and the second on temperature. 
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